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Week #4 Finy 
Modeling in the Time Domain........... eere nenne EE302 


As systems become more complex, representing them with differential equations or transfer 
functions becomes cumbersome. This is even truer if the system has multiple inputs and 
outputs. This lecture introduces the state space method which largely alleviates this problem. 


In contro! engineering, a state space representation is a mathematical model of a physical 
system as a set of input, output and state variables related by first-order differential equations 


The state space representation of a system replaces an nth order differential equation with a 
single first order matrix differential equation. For any system such as electrical or dynamical 
system, there exist a set of variables (states) such that if the values of the states are known an 


any reference time tOand if the input is known for t >t0then the output and the states can be 
determined for t =t0. 
LG 


State variable model 





1- The State Equations 
In the standard form the mathematical description of the system 1s expressed as a set of n 


coupled first-order ordinary differential equations, known as the state equations, in which the 
time derivative of each state variable is expressed in terms of the state variables 
x1(t),..., xn(t) and the system inputs ult), «in; n5). 


In the general case the form of the n state equations is: 


a, = fi(x,u) WT 
» ww 


E m f (x,u) S 


Ya = ว ั (10) 


For an LTI system of order n, and with r inputs, Eqs. (1) become a set of n coupled first-order 
linear differential equations with constant coefficients: 


ディ = 05324 + ut? + + in’, + Oty + ... + Du. 
ช่ ว = แผ ว 1 + 02332 + +. + (ว ท ซิ ท + imita Boss E Basta 
Ca E Puta + dete L i.: + Onmin + Onit + ... + Deu 


where the cozfficients aij and bij are constants that describe the system. 
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which may be summarized as: 
x= ค + Bu 


where the state vector x is a column vector of length n, the input vector u is a column vector 


of length r, A isan n X n square matrix of the constant coefficients aij, and B isan n X r 
matrix of the coefficients bij that weight the inputs. 


2- Output Equations 
An arbitrary output variable in a system of order n with r inputs may be written: 
y(t) = quay + 0525 4 Carla + dit +... + ds lr 


where the ci and di are constants. If a total of m system variables are defined as outputs, 
the m such equations tay be written as: 


31 = (yam + Mola + --- + Cinn + day, + ... + dae, 

Ya = Cati | C222 | ews | Canin | dziur | TT | doy Uy 
: E Y 

Un = Guía + Gots + d . + Cmnfn + mila coo... + dau. 


ーー 


y = Gx + Du 


where y is a column vector of the output variables yi(t), C is an m X n matrix of the constant 
coefficients cij that weight the state variables, and D is anm X r matrix of the constant 
coefficients dij that weight the system inputs. For many physical systems the matrix D is 


the null matrix, and the output equation reduces to a simple weighted combination of the 
state variables: 
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Summary: A system is represented in state space by the following equations: 


Input vector u Output vector y 





System 


u(t) Y 
| a descibed by state variables |...» i 









uc (t) 0x, ... x, } á y, の 
for t> to and initial conditions, x(t0), where 
X — state vector (nx 1) 


x = derivative of the state vector with respect to time (n x 1) 


y = output vector m f. E (m x i) 
u — input or control vector (r x 1) 
A — system matrix (n X n) 
B = input matrix | (n X r) 
C — output matrix (m xn) 
D — feedforward matrix | mx r) 


State variables: the smallest set of linearly independent system variables such that the values 
of the set members at time t along with known forcing functions completely determine the 
value of all system variables tor all t=t0.- 
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State vector: vector whose elements are the state variables. 
State space: The n-dimensional space whose bases are the state variables 


State equations: set of n simultaneous, first order differential equations with n variables 
where the n variables to be solved are the state variables 


Output equation: algebraic expression of the output variables of a system as linear 
combinations of the state variables and'the inputs. 


Selecting the state variables must follow these rules: 


1- A minimum number of state variables must be selected. 


The minimum number of state variables required equals the order of the differential 
equation describing the system. 





From the transfer: The order of the differential equation is the order of the denominator of 
the TF after canceling common factors in the numerator and denominator. 


A practical way to determine the number of state variables is to count the number of 
independent energy-storage elements in the system 


(capacitors and inductors in electrical syst.; masses and springs in mechanical Syst.). 
2- 'They must be linearly independent. 


Q) Wnat is meant by linear independence ? 


No variable in the set can be written as a linear sum of the other variables in the set. 


Recommended procedure: 
- Write the simple derivative eq for each energy-storage element. 


- Solve for each derivative term as a linear combination of system variables and the input. 
- Each differentiated variable is selected as a state variable. 


input. 


- All other system variables and the output are written in terms of the state variables and the 
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Ex: Find a state-space representation of the system shown in the figure if the output is the 
` current through the resistor. 





| i (t) idy ป 
v0(+) i) |R ^T C 


ー — 


Step 1 Label all the branch currents in the network: i,,i,,i, 


Step 2 Select the state variables: write the derivative equation for all energy-storage elements 
(L and C) and choose the differentiated quantities as the state variables (Vo 1) 





systems. 


Step 3 Write the right-had sides of “step 2” and the output equations as linear combinations of 
the state variables and the input. Network theory (Kirchhoff, etc.) will be used for the electric 






dv. . , | 
oe =i, = fio ip v(t) 





L9. 


di V; = fo(ve, โน w(t) 






State eq 
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Ex: Find the state equati 
ee equations for the electrical network in figure below 
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EX Write the state equation for the network 










4 Q Ya 2 H 





“ 
© Define the state variables as current through the inductor and voltage across the 


capacitors. Write two node equations containing capacitors and a loop equation 
containing the inductor. The state variables arev,,,v. . and 7 . 


Node equations are 











VI a a Le, + と し 
dy เ * P 
cl 「 + cl E dms RO a 
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ร 7, 4 า น์ ๒ -” ฟะ ร 29 Ya ZE 
aud the ioop equation is 
di : 
L a _ 
— 「 c2 cV d = 0 c Ir = 0.5v , = POP a 
or 


V. 
X;|-| 0 -2 2lJlv, |+l0 2M 


l 
i| [os 05 oz | lo ol? m 
y = es 2 inf. 
f Uc, 3- Side 
ง = E L J| t 1. ๓ น เปิ น 41. 
L 
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Ex: Find the state e i 
| quations for the translational mechanical System in figure below 





de _ 
dt 7 
dv, _ K VE K 
dt MM, 3° MM, 
dix, dx; dx 
Migg + Du + Kx, Ex =0 = 4v3 
=. dvz ー K K . 1 
-Kx + M2— a+ Kx; = f(t) dt Le” ไฟ * tm O 
In vector-matrix form, 
Xi 0 1 0 0 x) 0 
V1 u —K/M, —D/M; K/M, 0 Vi 0 
X2 i 0 0 0 1 x2 ¥ 0 Hu) 
úz j K/M 0 -K/M 0 v2 1/5; 
oV 2. 
Z|- Xi ——>5 Zi = = La 
๑ * Pp K 
Lua X = Le = ニー セー を t CR 
i ; ! 
27 e. e La = X, = du 
TE る 
Mr He 2 
LC e \ " . デュ ย 
ー E .p E O da. 
๑ = と » | Z3 
Za - La Ily 
. : K の L 
du E O ^d 
Ha i 
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EX: Represent the following transfer function in state space. Give your answer In và 
matrix form. 


] N-s/r EN J N-shn 





Writing the equations of motion 
(s^ +s+1)X,(5) — sX,(S) = F(s) 
= sX (s) + (s? +5+1)X,(5) — X,(s)=0 
— X,(s)+(s? +5+4+1)X,(s) =0 


Taking the inverse Laplace transform and simplifying, 
Xy =x x 2 + f 


Defining state variables, Zi, 


£j — Apr Lo = Xj, Z3 = X5, 24 5 Azor Zs = โร จ» Zg = X4 
Writing the state equations using the definition of the state variables and the 


inverse transform of the differential equation, 


e. 


<。 X= AA) AFA ES ZA ZA ES 


£3 = X) F £4 





อ LE จ . 
£, — X) Xj7X4— X4 b X4 = 2) 24 ZZ, 
Zs = X3 = Ze 
๑ อ ง E 
0 1 O 0 ๐ a 0 
Al th Y € 1 
. lo o o 10.0 0 
mms. gab] จ ง อ ย 
.r5. Wow bb 1 0 
ü 8-4 D eb si 0 
»=[0 © 0 O 1 0]= 
- emmees re . a - rg aa เน ร ーーーーーー ーーーーーー 
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ation to State Space 







Converting from the Differential Equ 
: ‘no differential equation, 
Ex Obtain the state equation in phase variable form for the following 1 q 


PY, d'y ladys Su(t) 


dr de E e three state variables as follows 
The differential equation 15 third order, thus there are three state Y: 


y= = y, A) = = y, ly = = y 


and the derivatives are 


Or in matrix form 

x 0 1 Olix | 10 
x, {=| 0 0 | 1 |x, (|0 เ ว (7) 
* | |-4 -3 -2| ๒ | [1 
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Extract the state Space using simulink] 





1- Save the differential equation in the following model 


ーーーーmーーーーーーーーーーーーーーーーーーーー ーーーー ーーーーーー 
ーーーーーーーーー ツ ーーー Ae 





EIDE Model | ^. ーー Sa 
File Edit it View | ‘Simulation ba Tod Help 





Bias) 
S&u8|X*Gesdqioc| epa fo ๆ เซ พ ๒ ต ๑ ไต 











y 
ln 4 
EP 
Quti 
| 
Ready เก ศร ชะ ห ท เศ ด คน แน ธร ร 0] 


2- in the command windows type the following command 


>> [A,B,C;D]-liiüaud( 5E, model) ก 2 
A= 

0 1 0 

0 0 1 

-4 -3 -2 
13 = 

0 

0 

1 
C= 

1 0 0 
1) = 

0 
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. i ate Space 
Converting from a Transfer Function to State 5] 


: : : 'rm in numerator 
1- Convertino a transfer function with constant term in 





uh li 
sann -order linear 
This case corresponds to a linear system that can be represented a 


differential equation with constant coefficients like: 








n-l j 
£N. x ay) +a,y(t) = bju(t) 
dt^" dt 


The classical Transfer Function (TF) representation of that system is obtained applying the 
Laplace Transform to the diff. eq. 


YG) _ 
U(s) 















b 


0 
n i. =} 4 4. 
$ E, “MS + Ao 





ป (ร ) Y(s) 


— G(s) 








e líthe State Space (SS) representation has to be obtained, a convenient way to select state 


variables is to choose the output, y(t), and its n-1 deiivatives as the state variables. This is 
called the phase-variable choice. 
A SS representation using the phase-v 


ariable choice for the state variables is said to be in the 
Controller Canonical Form (CCF). 


Differentiating] . 
both sides 
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2- Converting a transfer function with polvnom 





merator polynomial 
For example, consider a third order SISO system with second-order nu poly 









b,s? +b,5 +b, 


G(s)= 
(s) s) +a, ast d, 


First, separate the transfer function into two cascaded modules 
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We can write: 







qo 


| In the 
Ws ) = 


time domain 
--- at 





l 
3 2 | 
S "Pans" HOES, 





U(s) 






Ps) - Vos! +b,5+b,) W(s) 








*3 PE — 47 X - AX, — dX, + u 


X3 


v(t) = ใ พ ั 9+ bw(t)+ b w(t) 
= b,x, t bx, +b,x, 
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EX: Represent the following transfer function in state space. Give your answer in vector- 
matrix form. 






(s^ +35 +7) 
(s + 1)(s? + 5s + 4) 






The transfer function can be represented as a block diagram as follows: 


アリ 
— FERAM ———— 


ว | 
$465 49544 s 357 


Writing the differential equation for the first box: 


x4 6x4 9 x4 Ax =r(1) 





Defining the state variables: 


X — X P 
£ We 
peg “ Dros 
Thus Xq = X 
Y = X ZEN 
Xa = “ว ー lock MAR ra m 
x, =—4x —9 x— 6 x+ r(t) =-4x, —9x, —6x, + r(1) 
y=x+3x+7x=7x +3, +, 
From the second box, 
In vector-matrix form: 
0 1 0 Q 
x=|0 0 1 ม + 070 
[-b.02 259 1| 
pzi-83 1| X 
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MATLAB Code 


= Transfer function: 
den=conv([ 1 1],[ 154)? ee 
G=tf (num, den) nw Bau FOS + 4 


[A,B,C,D]=tf2ss (num,den) 


' convert to canonical form' 


Af-flipud (A); ล = 
Af=fliplr(Af) 
Bf-flipud (B) Er wold —4 
Cf-fliplr(C) 1 O 0 
1 O 
'statespace' 
Statespace = ss (Af,Bf,Cf,D) 
C = D 
1 3 7 
convert to canonical form Statespace 
Af = 
a = 
1 O Xl x2 <3 
1 O O x1 ๆ i り 
i dd -^ x2 0 0 1 
x3 -4 -9 -5 
Af = 
ไว == 
" " > x1 0 
-4 -9 —6 
xz O 
x3 1 
BE = 
ES = 
O 
Ù X1 x2 x3 
1 yl És 3 1 
cl = 
Cf = dt 
? 3 1 vi 0 






~ es 
— — — 
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Converting from a State Space t 


Given the state and output equations 
x = Ax + Bu 
y = Cx Du 
take the Laplace transform assuming Zero initia] conditions: 


sX(s) = AX(s) + BU(s) 
Y (s) = CX(s) + DU(s) 


Solving for X(s) in Eq. 
(sl — A)X(s) = BU(s) 
or 
X(s) = (s1— A) 'BU(s) 


where I is the identity matrix. 
Substituting Eq. ( ` .) into Eq. ¢ ~ t) yields 


y(s) = C(sI - A)’ BU(s) + DU(s) = [Clsk — A) B + DJU(s) 


We call the matrix [C(sl — A)! B + D] the transfer function matrix, since it relates 


— C(sI—A) ! B+D 
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Ex Convert the state space and Output equations to a transfer function 





State space model 


X (0 1 ic 0 
X; -2 -3 A + ] | 
เอ 0 


A 
X 
y- [1 0] 4 [0] u 
k D 
. Using the expression for derived transfer function 
Y(s) | ] ] 


U(s) 28 43542 7 (s 4- 1)(s +2) 
T(s)=C(sI-AJ B+D 


td IIT 
tl e 


| ฯ | 1 1 
= 9| [-5—— 
S -3542 S| «8 bae+2 





of 
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Ex: Convert the state space and output equatio 


4 —1.5 2 
c= | [X + u(i) 
4 () 
y = es 0.625 |x 


The state equation is converted to a transfer function using 
G(s)- CsI- A) B (1) 


where 


-4 -1.9 2 
A= _B=| | and C=[1.5 0.625] 
4 0 0 
Evaluating (sI— A) yields 
i s+4 1.5 
(sI-A)= 
-4 S 


Taking the inverse we obtain 


P 1 s. —L3 
E-A) q > 
si+4s+6|4 s+4 
Substituting all expressions into Eq. (1) yields 


35+5 


G(s SU M ae 
ts) so +4546 


MATLAB code 


2=[ -4 -1.5 : 4 0 ]; Transfer function: 
B=[ 2 ; 0 ]; ง ธ + 5 

๐ =[1.5 0.625]; | | | | |  ------------- 
D=0; ธา “2 + 4 ๒ +6 


[num,den]= ss2tf(A,B,C,D) ; 
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Home Work 3 ----...........-... แ ห 962 


(Find the state-space representation of the following systems (P (๐ ๆ 0 ๓ ๕ ss h TE 





US, 
case 1 10 12 12 เย 
cepe pm * 
vil) =) ไส = IHE IF ra 
Case 2 
Case 3 


— fix = 1 Nest = 
E ET Mysik 
3 ps dmm 37 UKE firs 1 N-s/m 
i = 1 N-s/m = ] N- 
xi [----- M: =2k8 i My = I kg NN 
UIUC 


K=2N/m 
— RNA 


E 本 レー なれ に 
_, Frictionless 
SW) 


- ar oo 


dalt, 
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~ 
vase 4 
bi 
(ase 5 
K¡=4 Nim um 
A 0 — (600 に = 
fy = 3N-shm| Mi =1 kg Ir 3 Im | M, = 2kg] fp, = 2 NesimE 
: | } 
Case 6 
NS) ] N-m-s/rad Tu) ex 
8 N-m-s/rad 
HE de Ee 
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